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ABSTRACT 

Quasars that allow the study of IGM He II are very rare, since they must be at high redshift along 
sightlines free of substantial hydrogen absorption, but recent work has dramatically expanded the 
number of such quasars known. We analyze two dozen higher-redshift (z — 3.1-3.9) low-resolution 
He II quasar spectra from HST and find that their He II Gunn-Peterson troughs suggest exclusion of 
very early and very late reionization models, favoring a reionization redshift of z ~ 3. Although the 
data quality is not sufficient to reveal details such as the expected redshift evolution of helium opacity, 
we obtain the first ensemble measure of helium opacity at high redshift averaged over many sightlines: 
t = 4.90 at z ~ 3.3. We also find that it would be very difficult to observe the IGM red wing of 
absorption from the beginning of He II reionization, but depending on the redshift of reionization and 
the size of ionization zones, it might be possible to do so in some objects with the current generation 
of UV spectrographs. 

Subject headings: galaxies: active — intergalactic medium — quasars: general — ultraviolet: galaxies 



1. INTRODUCTION 

Reionization of the intergalactic medium (IGM) is a 
major cosmological milestone in the evolution of the uni- 
verse. Although reioni zation of hydrog e n and He I oc- 
curre d at z > 6 (e.g., iFan et ail 120061 : iKomatsu et al.l 
l2010t ). reionization of He II was delayed. Stars, the ini- 
tial source of much of the ionizing radiation, do not pro- 
duce enough high-energy UV photons to ionize He II, 
and thus this process occurred only when there was 
substantial hard ionizing radiation from quasars. Both 
theoretically and observationally, He II reionization is 
thoug ht to have occurred gradually at redshift s ^3-4 
(e.g ISokasian et aLll2002t lAgafonova et al.|[2005l) . Even 
under the hard photoionizing conditions at z ~ 2-4, He II 
outnumbers H I by a factor of r\ ~ 50-100 (jZheng et al.1 
I2004bt iFechner et all 12006ft . and thus He II has much 
stronger opacity than hydrogen. As a result, the most 
sensitive and direct probe of this reionization era is the 
He II Gunn-Peterson effect, analogous to the effect in 
hydrogen where a substantially neutral IGM causes Lya 
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line blanketing to crea te a trough in quasar spectra 
(G unn fc PetersonllT965h . 

There is a variety of indirect evidence that He II 
reionization occurred near z ~ 3. Much of this 
work has tried to exploit the IGM heating associ- 
ated with this epoch, and its effect on the H I Lya 
forest, either by examining the increased temperature 
using various statistics r elated to therm a l broadening 
(e^. iSchave et al.l [20001: iLidz et al.l [20091 iBeckerlFaLl 
l201Gt ) , or the opacity changes associated with heatin g 
(e.g. JBernardi et al.lf200l iFaucher-Giguere et alJl200"l . 
While suggestive of reionization at z ~ 3, such studies are 
fraught wit h difficult systematics. contradicting ob serva- 
tions (e.g., iMcDonald et al.ll2001t iKim et alj|2002ft and 
difficult theoretical i nterpretations ()Bolton et aJj|2009t 
iMcQuinn et al.l I2009D . Another approach has been to 
use ion ratios to determine the ionizing UV background 
on either side of the He II Lyman limit (228 A), which 
sho uld change as th e helium opacity changes (although 
see iFurlanetta 120091 ) . Such measurements are also dif- 
ficult and suffer from substantial uncertainty, but they 
generally agree with a z ~ 3 reionization (e.g.. ISongailal 
\199§. lAgafonova et all 120071k 

Because of the large uncertainty associated with these 
indirect measurements, it is valuable to pursue direct 
measurements of the He II opacity, as an independent 
method of constraining the redshift of reionization. Un- 
like the hydroge n Gunn-Peterson trough, which saturates 
at sh i ~ 10" 5 (|Fan et all 120061 . although see iMesingerl 
2010), the He II Lya Gunn-Peterson trough can re- 
main sensitive to ion fractions ^Heii <^ 0.01, due to the 
later redshift of its reionization, the lower abundance of 
helium versus hydrogen, the shorter wavelength of the 
line (rpp ~ At, v o ,), and density fluctuations in the IGM 
(|McQuinnl 2009). In addition, the large fluctuations in 
the ionizing background from the rare nature of the lu- 
minous quasars thought to be responsible for He II reion- 
ization make transmission possible even during the early 
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stages of reionization (jFurlanettol l2009') . 

With He II Lya at 303.78 A, observations of the 
Gunn- Peterson trough at the relevant redshifts must take 
place in the far UV, and are thus very sensitive to low- 
redshift hydrogen absorption. Since observations in this 
regime must take place in space, much of this work has 
been done with the Hubble Space Telescope (HST), al- 
though observations of a few brighter targets were pos- 
sible with the Hopkins Ultraviolet Telescope (HUT) and 
the Far Ultraviolet Spectroscopic Explorer (FUSE). The 
recently installed Cosmic Origins Spectrograph (COS) 
on HST is proving especially well suited for this work. 
Due to intervening IGM hydrogen, only a few percent 
of randomly selected z ~ 3 quasar sightlincs arc clean 
(no substantial H I opac ity) down to the He II break 
l|M0ller fc Jakobsenlll99a ). Verified He II quasars are 
therefore rare, although the advent of the Galaxy Evo- 
lution Explorer (GALEX) UV maps covering most of 
the sky has recently gre atly increased selection efficiency 
([Syphers et al.ll2009alibl . hereafter S09a,b). The difficulty 
of probing the redshift regime interesting for He II reion- 
ization is increased because the number density of lu- 
minous quasars may decline exponentially in z = 3-4.5 
(I0smerlll982t ISchmidt et al.lfl995l : [Richards et all [20061: 
iBrusa et al.ll2009D . 

To date there are 27 quasars observed to be free of 
substantial absorption down to the He II Lya break 
(with a few others that have observable flux at the 
break but also substantial intervening absorption). Six 
of these have relatively high S/N spectral observations: 
Q030 2-003 (z = 3.29, iHogan et all Il997i iHeap et al l 
20001) . HE 2347-4342 (z = 2 90. ISmette et all 12001 



Zheng et al. | I2004bt IShull et al.l [2010D, HS 1700+6416 
(z = 2.72.lFechner et al.l l2006D. PKS1935-692 (z = 3.18, 
lAnderson et al.lll999D . SPSS J2346-0 016 (z = 3.5) and 
SDSS J1711+6052 (z = 3.8) Cboth IZheng"et~aII [20081 
hereafter Z08). Our team recently found 19 others, rang- 
ing from z ~ 3.1-3.8, with reconnaissance using the 
HST Advanced Camera for Surveys/Solar Blind Chan- 
nel (ACS/SBC) prism (S09a,b). The remaining two 
known He II quasa rs are HS 1157+3143 (z — 3.00, 
Reimers et all l2005h and SDSS J1614+4859 (z = 3.8, 



Zheng et al.ll2005lL the latter of which was reobserved in 



S09a at increased (though still not high) S/N. 

In section [21 we use this growing number of He II 
quasars to make spectral stacks in redshift bins, thereby 
averaging over sightline and individual object variations. 
This averaging is a vital step, since the sightline varia- 
tion is expected to be substantial during He II reioniza- 
tion due to the fact tha t the ionizing pho tons come from 
rare, luminous q uasars dFurla netto 2009) and the IGM is 
inhomogeneous (jMcQuinnl 12009) . Using these stacks, we 
can start to answer whether or not seemingly unusual 
absorption profiles, such as that of SDSS J1711+6052 
(Z08), are indicative of overall IGM evolution, or sim- 
ply happenstance intervening absorption along particu- 
lar sightlines. We discuss the feasibility of observing IGM 
red wing absorption, and comment briefly on apparent 
emission lines in our spectra. 

In section O we use our substantially expanded sam- 
ple of He II quasars from S09a,b to determine the IGM 
helium opacity. We consider fits to the spectral stacks of 
section [2] as well as aggregate statistics using pixel-by- 



pixel opacity measurements of individual objects. Under- 
standing the redshift evolution of He II sightlines would 
shed light on the progress of He II reionization, so we 
here begin this consideration. We find that while our 
current data can provide the first good estimate of the he- 
lium Gunn-Peterson opacity averaged over a broad high- 
rcdshift range, they give no conclusive evidence about 
any redshift evolution. The methodology we establish 
will be applicable to future higher S/N observations. We 
conclude in section [4] 

2. AVERAGED QUASAR SPECTRA AND THE PROXIMITY 

PROFILE 

We conducted surveys in HST cycles 15, 16, and 16 
supplemental (GO programs 10907, 11215, and 11982) 
searching for He II quasars with the ACS/SBC prism 
PR130L. This prism covers roughly 1250 A (R = 380) 
to 1850 A ( R = 40). The spect ra were extracted us- 
ing aXe 1.6 (|Kiimmel et al.l l2009f ): further details may 
be found in S09a,b. The targeted quasars were identi- 
fied in the Sloan Digital Sky Survey (SDSS; lYork et al.l 
2000]), and detected as UV sources in the GALEX surveys 
(the all-sky AIS, the med ium-depth MIS, or the deep 
PIS: iMorrissey et al.1 120071 ). The cross-correlation cata- 
log thus created is presented and described in S09a. In 
this section we consider methods of combining these spec- 
tra to extract information about typical He II quasars. 

2.1. Stacks of Spectra 

Our ~4 ks HST exposures were intended as reconnais- 
sance to verify flux all the way down to He II Lya, but 
stacked spectra and a few higher flux quasars allow more 
detailed analysis, which we pursue here. The individual 
spectra for all objects can be found in S09a,b, along with 
details on the observations and data reduction. We found 
29 UV-bright quasars at z > 3.1, at least 23 of which have 
flux breaks at He II Lya, and 19 of which are free of any 
observed absorption redward of the He II break. Of those 
whose flux breaks redward of He II Lya, one has a clear 
continuum break to zero detectable flux ~15 A (rest) 
redward of the expected 304 A break, likely due to a 
low-redshift hydrogen Lyman-limit system (LLS), while 
the others have more extended, complicated absorption, 
perhaps arising from hydrogen absorption line blanketing 
near the Lyman limit of an LLS or a damped Lya ab- 
sorber (DLA) . Two of the 23 with breaks have relatively 
strong emission lines, with weak, low-S/N continua. Z08 
contains two other high-S/N ACS prism spectra that we 
use in our analysis. 

In order to average over sightline variance and extract 
information from the aggregate spectra, we coadd the 
spectra in redshift bins. We performed an observed- 
frame stack of these spectra in S09b, to verify that 
ACS/SBC PR130L does not have any strong instrumen- 
tal features uncorrected in the reduction. In S09a, we 
presented median rest-frame stacks of the 31 far- UV- 
bright spectra in three large redshift bins (3.1 < z < 3.3, 
3.3 < z < 3.6, and 3.6 < z < 4.1), finding no obvious 
redshift evolution, but pursuing no detailed analysis. 

In this paper, we improve the coadding method, and 
present a more detailed analysis of the resultant stacks. 
The resolution of ACS spectra depends in a highly non- 
linear manner on wavelength. Because each quasar is 
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at a distinct redshift, the restframe wavelength associ- 
ated with each pixel (the knots of a spline fit), as well 
as the restframe wavelength range spanned by each such 
pixel (i.e., the dispersion solutions transformed to the 
restframes) , will differ for each quasar to be combined in 
the spectral stack. We must, therefore, map the spectra 
to a common wavelength grid before coadding. In S09a, 
we did this by spline fitting spectra to the wavelength 
knots of the highest redshift (lowest resolution) quasar. 
This is not a bad approximation, but is suboptimal for a 
number of reasons. It limits the upper wavelength in the 
rest frame to that of the highest- redshift object, ignoring 
data we have beyond that; it neglects any finer detail that 
we can extract from the higher-resolution spectra, since 
it uses the lowest-resolution wavelength knots; in some 
circumstances, there is not an obvious way to extract an 
error spectrum for the resultant stack; and finally, while 
the spline-fitting method is good for wavelength knots 
near each other, it is somewhat worse for the substan- 
tially shifted knots we have here. 

Our approach in the current work is to make the spec- 
tra almost everywhere continuous by replacing each dis- 
crete spectrum by the sum of a set of top-hat functions, 
one for each pixel. We then define a wavelength grid, and 
for each pixel of this grid, calculate the overlap with every 
pixel in the individual spectra (adding all contributions 
in each individual spectrum, weighted by their overlaps). 
We then take the median over all the input spectra for 
each pixel in the stacked spectrum. We considered other 
methods, such as taking the mean or a S/N- weighed 
mean instead of the median, but we found in simulations 
of noisy spectra that the median is marginally better 
even when there is no sightline variance. With sightlinc 
variance, the median is certainly more robust. 

Prior to this coadding, we shift the spectra to their 
respective rest frames, using the SDSS pipeline redshifts 
(verified by eye to be accurate, except in one case where 
we shifted an object by Az = 0.05). Then we normalize 
each spectrum by its continuum flux (avoiding any broad 
rise towards He II Lya) . In the rest frame, we stack the 
25 spectra with He II Lya breaks, using all ACS spectra 
in S09a,b and Z08 with significant flux breaks, except 
those with LLS (since they have extremely low S/N con- 
tinua, and are heavily affected by hydrogen absorption) . 
We separate the spectra into four evenly-spaced redshift 
bins: 3.1 < z < 3.3 (7 objects), 3.3 < z < 3.5 (7 objects), 
3.5 < z < 3.7 (5 objects), 3.7 < z < 3.9 (6 objects). Not 
all spectra contribute to every wavelength, but all con- 
tribute at the He II Lya break. To find the error on the 
median, we use a bootstrap method over the medians 
of all contributing spectra at each wavelength, smoothed 
with a running average over three pixels; for further dis- 
cussion of our bootstrap methods, see section 13.11 The 
stacks are shown in figure [TJ 

2.2. Fits of Spectra 

We fit a simple model quasar spectrum with an ab- 
sorption edge to the stacked spectra and a few higher 
S/N individual objects. We use a power law continuum, 
with Gaussian emission and absorption lines, and an ab- 
sorption edge. There is not expected to be any change 
in spectral index over our limited wavelength region, and 
due to high Doppler parameters, the emission lines are 
approximately Gaussian rather than Voigt profiles, so 
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Fig. 1. — Stacked spectra in redshifts bins with fits overplotted 
(dashed red lines and spectrum errors shown (dotted lines). There 
are 7 objects in redshift 3.1-3.3, 7 objects in redshift 3.3-3.5, 5 
objects in redshift 3.5—3.7, and 6 objects in redshift 3.7-3.9. 

this should provide a good basic fit. We model the IGM 
absorption as an optical depth r(z) = tq/(1 + w(z)), 
where To is the optical depth far from the quasar (where 
the intergalactic UV background dominates), and w is 
the ratio of ionizing flux from the quasar and from the 
intergalactic background. To determine t , one must 
have observations far from the quasar, at w <C 1, which 
might be possible with some (though not all) of our spec- 
tra. We take the simple model w{z\) = wo/(Dj j (zi, 22)), 
where is the luminosity distance from the quasar (in 
practice, from the absorption edge at Z2). Physically, 
wq depends on the luminosity and SED of the quasar, 
and thus is of scientific interest for future high-quality 
spectra. However, the current data do not allow us to 
robustly determine its value, and therefore we treat it 
as a nuisance parameter that we allow to vary in our 
fits, constraining it only by a quasar lifetime limit of 10 8 
years, w hich should be appro priate for the se luminous 
quasa rs ()Hopkins et al.l 120061 although see iKellv et all 
I2010T) . Simulations show that this geometric dilution 
factor should be a good fit to the proximit y zone ev- 
eryw here except < 1 Mpc from the quasar (jPartl et al.l 
2010), which is not resolvable in our spectra. We use 
a ACDM co smology with h = .704, f2 A = 0.73, and 
fi TO = 0.27 ()Komatsu et al.ll2010D . although our results 
are insensitive to these values, since we do treat Wq as a 
nuisance parameter. 

We fit this simple model to our spectra using the 
Leve nberg-Marquardt l east-squares fitting program MP- 
FIT (jMarkwardtl 120081 ). As a first step, we fit the ob- 
served data with just the analytic quasar spectrum de- 
scribed above. We use this to find the line centers, and 
the wavelength of the absorption edge, and obtain initial 
values for line widths and amplitudes. In the second step, 
we fit the data using convolution of the model through 
the ACS prism's (nonlinear) instrumental resolution. In- 
cluding this convolution makes the fitting step much less 
efficient, so we fix the wavelengths of line centers and 
the absorption edge at this point, as they are negligibly 
affected by instrumental resolution. 

The fits for the S09a,b targets 1006+3705 {z = 3.2) 
and 1253+6817 (z — 3.47) are shown in figures |2^,b. 
We also include fits for the Z08 targets 1711+6052 and 
2346-0016 in figures [Sfc, d. The first three fits use no ab- 
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Fig. 2. — Spectra of four individual quasars with fits overplotted 
(dashed lines, red in online version) and spectrum errors shown 
(dotted lines), (a) SDSS J1006+3705, z = 3.20, fit with six emis- 
sion lines, including the only use of He II Lya in any of our fits in 
this figure or figure[T] Reduced x 2 = 1.15. (b) SDSS J1253+68I7, 
z = 3.47, with with one emission line. Reduced x 2 = 0.97. (c) 
SDSS J1711+6052, z = 3.82, with two emission lines. Reduced 
X 2 = 0.89. (d) SDSS J2346-00I6, z = 3.51, with five absorption 
lines. Reduced x 2 = 0.68. 

sorption lines, while the 2346 fit uses only absorption 
lines, with no emission lines. Below we comment briefly 
on the apparent emission lines, and we consider the IGM 
opacities of these fits in section 13.11 In every case, the 
break wavelength found in our fit is consistent with He II 
at the quasar redshift. The redshift-binned stacks are fit 
in figure [TJ 

Due to the high opacity of the IGM at A < 912 A 
from H I continuum absorption, clear sightlines to ex- 
amine quasars in the extreme UV (EUV) are rare. Our 
own galaxy's opacity constrains these sightlines to be at 
higher redshift (z > 912 A/A - 1). In the process of 
finding He II quasars, we discovered a number of such 
clear sightlines appropriate for EUV quasar spectroscopy. 
Some quasars appear to have multiple observable emis- 
sion lines in this region, (see, for example, figure [2^i), 
although it is interesting that these lines may not be as 
universal in line strength as the well-known lines seen in 
the optical and the near and far UV. 

We have run simple broad emission line region mod- 
els using the photoionization sof tware Cloudy 08.00 (last 
described in lFerland et alJ fl998). with a s tandard range 
of de nsity and ionization paramters (e.g.. iKorista et all 
1997), so lar and supersolar metall i cites, and a variety 
of SEDs dMathews fc Fe"rlaiidlll987t IKorista et all 119971: 
IHaro-Corzo et al.1 120071) .~ These models indicate that 
the strongest emission line in this portion of the EUV 
(A ~ 304-430 A) is He II Lya, with an equivalent width 
of ~6-ll A for a covering factor of 10%. This strength 
is easily resolvable in our spectra, and the line is pre- 
dicted to dominate by one or two orders of magnitude 
over any other lines, so it is therefore puzzling that our 
averaged and individual quasar spectra (figures Q] and [2]) 
do not show obvious He II Lya, with the exception of 
1006+3705. 

Our spectra may show a noticeable emission line near 
313 A, where photoionization models predict a C IV line. 
Although this line is predicted to be stronger than most 
other EUV lines, the difference is fairly small, and the 
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Fig. 3. — Example of how an IGM red absorption wing affects a 
Gaussian He II Lya emission line, for a z = 3.8 quasar (assumed 
to be near the beginning of He II reionization, and thus having 
the strongest possible red wing). The red wing might substantially 
reduce the line's EW and peak flux relative to the case with ab- 
sorption only below the break (the EW decreases 26% with the 
inclusion of the red wing, in this example). However, this is obser- 
vationally degenerate with an intrinsically weaker line, unless the 
very narrow absorption profile can be seen (and distinguished from 
a proximity profile). The low-resolution FUV grating of COS has 
a dispersion of 0.56 A per resolution element; the highest-redshift 
He II quasars currently known are too faint for observation with 
COS's medium-resolution grating. 

models uniformly predict this line to be substantially 
weaker than He II Lya. We find that microturbulence 
is an effective (and plausible; e.g.. iBaldwin et afll2004l ) 
way of increasing metal line strength relative to He in 
the EUV (as has been previo usly noted for metals com- 
pared to H in the FUV/NUV: fFerland 11951 . but it still 
leaves the strongest metal lines over ten times weaker 
than He II Lya. The four well-known He II quasars 
that have been observed at higher spectral resolution 
and S/N have not obviously shown any metal lines, or 
even He II Lya consistently. Further analysis of poten- 
tial emission lines in our ACS spectra is ongoing, but 
it appears that definitive answers await the many He II 
quasar spectra our team and others are obtaining with 
the higher resolution of COS. 

We do see an absorption profile reminiscent of the fea- 
ture in 1711+6052 (figure [5b) in some of our individual 
spectra (e.g., figure!!]}) and many of our averaged stacks, 
albeit generally not one as prominent as that of 1711. 
One possible explanation is that He II Lya is absorbed 
by the IGM, since our quasars are chosen to lie in an 
era with a substantial He II fraction. A neutral (or for 
helium, singly-ionized) IGM is expected to produce a red 
wing of absorption near resonant lines ([Miralda-Escudel 
1998), as in figure This is a possibility of great inter- 
est, as we now have several known high-redshift (z > 3.7) 
He II quasars (Z08, S09a,b), which may exist in an IGM 
with a substantial He II fraction. For XHeii ~ 1, the 
Gunn- Peterson t rough would be saturated (although see 
lFurlanettoll2009[) , but the red damping wing is most eas- 
ily observed in just these conditions, and would be a way 
of directly constraining high x hrtt- For our red win g pre- 
dictions, we take Of, = 0.0456 (iKomatsu et al.ll2010f) and 
helium mass fraction Y p = 0.2483 (Stcigman 2003). 

However, as one might expect, we predict that such a 
red wing would be extremely difficult to detect in He II. 
When observing an object with an ionized proximity 
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zone, such as many quasars, the red wing is expected 
to be subst antially weakened and thus difficult to detec t 
in H I Lya (|Madau fc Reedl2000t iCen fc Hainrnn1l2000l ). 
For He II Lya, however, this red wing should be reduced 
even compared to H I Lya, largely because of the lower 
abundance of helium and the later redshift of its reion- 
ization (leading to a lower density). 

A possible factor intensifying the red wing is that 
quasars probably form in denser regions, while analytical 
red-wing analyses have assumed a uniform IGM. While 
moderate overdensities might still be insufficient to pro- 
duce an observable red wing, substantial overdensities 
(such as He II DLAs) might be strong enough to do so 
(for an example possibly involving a high-redshift galaxy, 
see Z08). Denser environments would have higher ab- 
sorption because of the higher density, but they also have 
more photon sources, and thus might be more highly ion- 
ized, reducing absorption. Apart from the quasar prox- 
imity zone itself, however, this increase in ionization is 
probably not significant, because quasars are by far the 
dominant producers of He II-ionizing photons. 

While any diffuse IGM He II red wing might be dif- 
ficult to detect in an absorption profile (and impossible 
with the resolution of our present ACS spectra), it will 
likely reduce the flux of the He II Lya line (see figure 
[3]). Of course, this effect could not be discerned unless 
the intrinsic line flux were known, but it does reduce 
the He II Lya strength relative to our photoionization 
models. However, this reduction is moderate, even with- 
out an ionization zone diminishing the red wing, and 
could not cause the line to disappear. The red wing effect 
should be larger at higher redshift, both because the IGM 
helium is less ionized and because the IGM density is 
larger. However, we see in SDSS J1137+6237 [z = 3.78; 
S09b) what certainly appears to be a He II Lya emission 
line that rises right until the He II Gunn-Peterson break. 
This confirms our calculation that IGM red wing effects, 
if present, are small and not resolved by the ACS/SBC 
prism. 

It might be possible to observe this red wing absorp- 
tion with COS, given the fortuitous circumstances of a 
relatively bright He II quasar at high enough redshift 
to be near the onset of He II reionization, and with- 
out a substantial proximity zone. Such requirements 
are not impossible to meet. For example, the very lu- 
minous quasar H E 2347-4342 shows no proximity zone 
(Shull et al.ll20ld ). although it is at too low a redshift to 
be of use for this test. The known high-redshift He II 
quasars are quite faint, and require using COS's low- 
resolution grating, which has a dispersion of ~ 0.6 A 
(barely able to distinguish the red wing from average- 
density IGM, in a noise- free spectrum; see figure [3]). 
Brighter quasars, should some be discovered at z > 3.7, 
or should the epoch of He II reionization start later 
than expected, could use the medium-resolution grating 
of COS, with a dispersion of 0.07 A. One advantage of 
looking for the red damping wing during the He II reion- 
ization epoch, compared to that of H I, is that there is a 
known substantial quasar population at the presumptive 
beginning of helium reionization (z ~ 4), but not at the 
beginning of hydrogen reionization. 

There is another IGM absorption effect reducing 
He II Lya emission that does not require a strong red 



wing. It is known that the high-ionization emission lines 
of quasars are syste matically blueshifted relative to low - 
ionization lines (e.g..lGaskcll 19821 iRichards et aLll2002ft . 
and thus presumably also relative to the systematic red- 
shift of the quasar an d neighboring IGM. For example, 
iRichards etall (|2002l ) find that C IV 1549 is shifted by 
a median value of nearly 1000 km s -1 relative to Mg II 
2799, and in certain instances the shift may be up to 
~3000 km s . Shifts such as this would move He II Lya 
emission line substantially into the region covered by 
IGM He II absorption, without relying on the red wing 
for absorption. However, the same ionized proximity 
zone that destroys the red wing also impacts this effect. 
It should be noted that of the four well-studied He II 
quasars, Q0302-003 and PKS1935-692 have clear p roxim- 
ity zones (Heap et al. 2000; Anderson et al.lll999D while 
HS 1700+6416 and HE 2347-4342 do not (iFechner et al.l 
120061: ISmette et al.l[200l IShull et alJf2010h . Despite this 
difference, none of the four have clear He II Lya emission 
lines, even in higher-resolution spectra. HS 1700+6416, 
at z — 2.72, is in an IGM where He II is substantially 
reionized, and thus the IGM opacity is much smaller. 
Even so, the quasar does not show clear He II Lya emis- 
sion. Thus, while the blueshifting of emission lines might 
be a contributing factor to the weakness of He II Lya in 
some objects, it seems unlikely to explain the whole effect 
alone. 

Because of the small sample size, there is always the 
possibility that some of the variation we see in emission 
lines and absorption profiles is just due to intrinsic dif- 
ferences between quasars. We do see some objects with 
apparent He II Lya emission, sometimes strong (e.g., 
SDSS J1442+0920 and SDSS J1137+6237, both S09b, 
and SDSS J1315+4856, S09a). We have enough objects 
that we know strong He II Lya emission is not com- 
mon, as the stacks of section 12.11 show, but we await 
higher-resolution COS spectra to distinguish the effects 
of emission and absorption. 

3. IGM HE II OPACITY 

In this section we consider some ways of combining 
our spectra to extract information about the IGM He II 
opacity. At the most basic level, our spectral data al- 
low for helium opacity measures along several individual 
sightlines, as well as an ensemble measure of the me- 
dian optical depth averaged over many sightlines. With 
a large number of high-redshift sightlines clean down to 
He II, it is also natural to see if our data reveal any red- 
shift evolution of the He II Gunn-Peterson opacity. There 
is some evidence from H I Lya forest observations and 
inferred temperatures that He II reionization occurred 
earlier than some previous studies suggested, at z ~ 3.4 
dLidz et al.l l2009( ). Although this result is interesting, 
Lidz et al. (2009) caution that it is statistically marginal, 
and so it would be useful to directly examine the Gunn- 
Peterson opacity to constrain this. The ACS/SBC spec- 
tra taken in our reconnaissance programs were intended 
to verify clean sightlines only down to He II Lya, and as 
a result, some of our objects have very little data below 
the He II break. While detailed analysis of the Gunn- 
Peterson trough in those particular objects is not pos- 
sible, our large sample includes many objects at higher 
redshift that do allow for such an analysis. 
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3.1. Hell Gunn-Peterson Trough 

The first measure of opacity we turn to is that obtained 
in our spectrum fits of individual objects (figure [2]) . No 
region of our spectrum for 1006+3705 (z — 3.20) is in 
the w <C 1 case (i.e., where the IGM ionizing background 
dominates over the quasar flux), and we therefore have 
no good determination of To. There is somewhat more 
hope for 1253+6817 (z — 3.47), as we have a long trough 
(although with low S/N). Our fit finds To ~ 5 far from 
this quasar, at a mean redshift of z ~ 3.25, while at a 
similar redshift, 2346-0016 (z = 3.51) is fit with to ~ 6. 
We fit 1711+6052 (z = 3.82) with t ~ 4, lower than 
those at lower redshift, but this difference is unlikely to 
be significant. There are large systematic uncertainties 
in these measurements, arising both from our choice of 
model and from the expected substa ntial sightline vari- 
ance (e.g., iFurlanetto fc Dixonll2010l ). and we therefore 
present these values only as a description of our fits, and 
not as reliable measurements of the IGM opacity. 

Fitting the stacks (figure [1]) provides another estimate 
of the IGM opacity, one that averages over sightline vari- 
ance. We formally obtain fits for the opacity that are 
not unreasonable (to ~ 4-6, consistent with our more 
robust estimates found below) , but the systematic errors 
associated with these determinations are still very large, 
and limit their usefulness. The model we are fitting to 
the data is also simplistic, in that it assumes a constant 
opacity in the Gunn-Peterson trough of any given spec- 
trum, with no redshift evolution. We turn, therefore, to 
our entire data set as a collection of many opacity data 
points, allowing for redshift evolution, and averaging over 
individual objects, spectrum extractions, and exposures. 

We pursue here an analysis of average He II Gunn- 
Peterson opacity in the redshift range z = 3.0-3.7, 
binned by Az = .1 to look for redshift evolution. 
iDixon fc Furlanettol (|2009T ) looked for evolution in the 
range z = 2.0-3.2, using literature values for the hand- 
ful of well-known sightlines with data covering some of 
this region (at most redshifts in this range, only one or 
two quasars contributed data). While we have the ad- 
vantage of a far greater number of sightlines, we also 
have the disadvantage of much lower-resolution, lower- 
S/N spectra, although the inclusion of the longer expo- 
sures of Z08 helps somewhat with the S/N. (We note, 
however, that an analysis done without these two long 
exposures still gives results consistent with those pre- 
sented here.) The background subtraction for the ACS 
prism, while generally acceptable, is a source of consider- 
able uncertainty when extracting such a small signal. In 
creating their lower - redsh ift opacity evolution analysis, 
IDixon fc Furlanettol (|2009| ) were able to ignore those lit- 
erature opacity values that were only lower limits. How- 
ever, due to the relatively high noise of our short expo- 
sures, as well as the higher opacity one would expect at 
this higher redshift regime (prior to the completion of 
He II reionization), we cannot ignore data points that 
only supply lower bounds on the opacity. 

For each object, we define a portion of the spectrum 
to be used as its Gunn-Peterson trough, excluding the 
very noisy data at A < 1250 A, as well as any obvious 
ionization zones near the quasar's He II Lya break. We 
fit a power law to the continuum region (avoiding the 
region around any He II Lya emission), and extrapolate 



this to find the expected continuum flux in the absence 
of IGM helium absorption, / C (A). Each data point in the 
Gunn-Peterson trough has a signal flux, /gp (after the 
background is subtracted during spectrum extraction). 
We find the ratio /gp(A)// c (A), and average all ratios in 
each redshift bin. By working with the ratio, we avoid 
absolute flux calibration problems, which might be sig- 
nificant for faint objects in these prism exposures (S09a). 
This approach also avoids the highly skewed distribution 
of the optical depth, where, because of large or effec- 
tively infinite opacities, the mean and median are very 
different, with the former generally not defined. Even 
the ratio distribution is skewed, however, so we prefer 
the median as a much more robust estimator than the 
mean. 

Because the signal ranges from being of the order of 
the background to much smaller, we have a number of 
unphysical data points for which /gp < 0, simply due to 
fluctuations in the background. In order to not bias the 
average, a standard procedure is to include these unphys- 
ical values while taking the average fe.g.. lCowan1ll998L 
pp. 136-142), and we use this for our point estimates of 
the opacity. However, with lower numbers of data points, 
and /gp ~ 0, even the average can attain unphysical val- 
ues. In this case, no point estimate exists for the opac- 
ity, but we would still like to obtain one-sided confidence 
intervals. The question of how to find reasonable confi- 
dence intervals for data near physical boundaries is quite 
difficult, even for very simple probability density func- 
tions (p.d.f.s) (e.g., for a review of methods to deal with 
Gaussian and Poissonian data nea r boundaries, and their 
difficulties, see Man delkernll200l . 

Since the p.d.f. of our opacity is both complicated and 
substantially unknown, we take a non-parametric boot- 
strap Monte Carlo approach. To ensure adequate data 
for resampling, we require n > 20 data points in each red- 
shift bin (in practice we have n > 24), which means that 
we must discard the poorly populated z = 3.7-3.8 bin. 
In each bin we perform N = 10 5 bootstrap realizations. 
Because the log function is monotonic, the confidence 
interval found for the ratio distribution transforms with 
equivalent cover to an interval for opacity. Note, how- 
ever, that while this transformation invariance is true 
for standard bootstrap percentile con fidence intervals, it 
is no t true for percentile^ methods (|DiCiccio fc Efronl 
119961 ). Thus despite the latter technique being better 
in some circumstances, we avoid it here. The statistical 
analysis we perform has some similariti es to that inde- 
pendently created by iPress et al.l (|1993l) . although they 
were not working near the zero bound and could discard 
negative flux measurements. 

In standard bootstrapping, which is a technique of re- 
sampling with replacement, every data set generated by 
the bootstrap contains n data points, where n is the num- 
ber of original data points observed. The statistics of 
this ensemble of mock data sets can then be calculated, 
and assumed to approximate the unknown p.d.f. A sim- 
ple but quite broadly applicable method of protecting 
the bootstrap against failure is the m-out-of-n method 
of bootstrapping, where each boo tstrap data set is gen- 
erated with m < n data points (|Chernic a i2008h . The 
required asymptotic behavior of m(n) is that as n — >• oo, 
Tti — y oo and m/n — > (e.g., m = n a , for < a < 1). In 
practice, our specific choice of m is found via the mini- 
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mum volatility method (jPolitis et al.lll999l pp. 197-200). 
While our data are not subject to most of the problems 
the m-out-of-n method addresses, there is at least one 
case in which it might be helpful. For a signal at zero, 
one expects about half of the data points to lie in the 
unphysical negative regime (after background subtrac- 
tion). Although the median is a fairly robust estimator, 
it has a breakdown point of 50%, and so when it happens 
that more than half the points in a bin are unphysical, 
the median itself becomes unphysical. This problem oc- 
curs in our z = 3.1-3.2 bin, where ^53% of the data 
are unphysical — not enough to think that this is due to 
anything other than a random fluctuation (18% proba- 
bility) around a very small (consistent with zero) signal, 
but enough to interfere with the determination of our 
confidence intervals. In this case, confidence intervals 
with larger coverage help, but so does the m-out-of-n 
approach. We conservatively display m-out-of-n confi- 
dence intervals for all redshift bins, although the differ- 
ence from the n-out-of-n intervals is insignificant in most 
cases. 

Figure |4] shows Gunn-Pctcrson opacity versus redshift, 
and it is clear that no statistically significant evolution 
can be seen. Where available, the median estimate is 
plotted as an 'x', and the confidence limits as arrows (off- 
set slightly in the horizontal direction for clarity). The 
thicker black arrows indicate standard bootstrap lower 
limits, 68% for those on the left and 95% for those on 
the right, while the thinner blue lines show 68% and 95% 
confidence limits for m-out-of-n bootstrap estimates. To 
enable comparison between redshift bins, only one-sided 
1 — a limits of the opacity are plotted, so in the cases 
where finite two-sided confidence intervals exist, when 
one combines the upper and lower limits shown it pro- 
duces a 1 — 2a two-sided confidence interval. While the 
z = 3.1-3.2 bin looks anomalous, this appears to be due 
to the data problems discussed above. The z = 3.2-3.3 
bin has upper limits defined for most confidence levels, 
and we can therefore state that when considering the 
two-sided 90% intervals, there is no significant difference 
between this bin and that of z = 3.1-3.2; i.e., the ap- 
parent anomaly is not statistically signficant. We con- 
clude that our data are unable to reveal any evidence 
of opacity evolution in this redshift regime (a conclusion 
confirmed by shifting bin edges and sizes), although we 
caution that this says more about our data than it does 
about the IGM — note that three of the six bins (includ- 
ing those at highest redshift) include infinite opacity in 
their two-sided 90% confidence intervals. 

We aggregate the data further to make an opacity esti- 
mate for this entire high redshift region, which is mostly 
inaccessible to the well-known He II quasar sightlines. In 
the large bin z = 3.1-3.7, we find a median r = 4.90+^4 
(95% confidence, m-out-of-n). Because the ACS prism 
provides more data points at lower wavelength (and thus 
redshift), the mean redshift of the data used for this 
opacity is z w 3.3. This opacity at this redshift agrees 
well with the r = 4.75j j^ measuremen t of Q0302-003, 
found for z w 3.1-3.2 (jHeap et al.l 120001 ). and the weak 
lower limits of r = 6^^° for z w 3.4-3.5 (jZheng et al.l 
I2004al) and r > 3.1 (90% confidence) for z 3.35 (Z08). 
Comparing our binned data of figure |H to the mo dels 
considered in figure 7 of iDixon fc Furlanettol (|2009l ). we 




Fig. 4. — He II Gunn-Peterson opacity as a function of redshift, 
as derived from our ACS spectra. Where available, median esti- 
mates are plotted as 'x's, and confidence interval limits are plotted 
slightly offset in the horizontal direction. In each bin, where avail- 
able, offset to the left are the 68% confidence limits, and offset to 
the right are the 95% limits. Upper and lower limits are calcu- 
lated independently; note that when combining independent 1 — a 
upper and lower limits, the resulting confidence interval is 1 — 2a. 
The thick arrowheads denote standard bootstrapping limits, while 
the thinner blue tails are for m-out-of-n methods (see text). The 
apparently anomalous opacity of the lowest redshift bin is not a 
statistically significant deviation from the other bins (see text for 
a discussion). 

find that our measured opacity is consistent with their 
model of He II reionization ending (full ionization) at 
^roi ~ 3, and not consistent with the extreme model 
of z ro i > 3.8 (it is also inconsistent with an extrapola- 
tion of the other extreme model, z ie i = 2.5). We cau- 
tion, however, that the models of IDixon fc Furlanettol 
(2009) require assumptions about the spectral index of 
quasars in the EUV. This is poorly k nown at A < 912 A 
(jTelfer et al.l 120021 ; iScott et al.ll2004D . and has no direct 
observational constraint at the energies needed to ionize 
He II, A < 228 A. This uncertainty relaxes the constraints 
we can place o n the redshift of full helium reionization 
(IMeiksinll2005l) . 

The ACS spectra here, with the exception of those from 
Z08, were obtained with the intent of verifying the exis- 
tence of He II Gunn-Peterson troughs. As such, it was 
not expected that they themselves would be a partic- 
ularly sensitive probe of the He II opacity, but it was 
worthwhile investigating this question, as they consti- 
tute by far the largest sample of He II quasars to date. 
Indeed, beyond z ~ 3.2 the only opacity data is from 
these ACS spectra (includin g the Z08 objects in o ur sam- 
ple) and the weak limit of iZheng et all (|2004a[) (which 
was from SDSS J2346-0016, reobserved and included at 
higher S /N in our sample) . It is therefore perhaps disap- 
pointing, but hardly surprising, that the uncertainties in 
the current data do not allow any robust determination 
of the He II opacity evolution in this interesting redshift 
regime. Higher S/N spectra of the brighter of our quasars 
will be invaluable for making this direct determination of 
opacity; a number of such observations are now planned 
for HST. 

4. CONCLUSION 

The relatively large numbers of He II quasars our team 
has recently found allow us to construct average spectra 
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in four high-redshift bins. While the individual objects 
are often low S/N, one of the primary sources of uncer- 
tainty in He II Gunn-Peterson studies is sightline vari- 
ance, which we can now begin to average over. We find 
that the absorption profile of SDSS J1711+6052 may not 
be as unusual as was initially thought. We cannot resolve 
the puzzle of why He II Lya emission appears so weak in 
many (though not all) known He II quasars, but we con- 
sider two possibilities, and find that neither a red wing of 
IGM absorption nor emission-line blueshifting can alone 
account for this. It is unlikely to find He II quasars that 
will allow observations of the diffuse IGM red wing at the 
beginning of He II reionization, but denser IGM knots 
and rare young quasars might provide such opportuni- 
ties. 

We find an IGM He II r GP - 5 at z = 3.1-3.7, in broad 
agreement with the few other probes of this high-redshift 
regime. Our data are not of sufficient quality to constrain 
the redshift evolution of this opacity to any significant 
degree, although we may be able to rule out models in 
which He II reionization ends extremely early (z > 3.8) 
or extremely late (z ~ 2.5), assuming particular quasar 



SEDs. Estimates of the Gunn-Peterson opacity derived 
from model fits of single, higher-S/N quasars agree with 
those found from ensemble statistics. 

The discovery of more He II quasars will allow fur- 
ther reduction in sample variance, and higher-resolution, 
higher-S/N observations with COS will allow better con- 
straint on both absorption profiles and He II opacities. 
Observations with both these aims are currently under- 
way. 
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